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Abstract 

By using the gauge/gravity duality and the Maldacena prescription we com- 
pute the expectation values of the Wilson loops in hot, noncommutative Yang-Mills 
(NCYM) theory in (3+1) dimensions. We consider both the time-like and the light- 
like Wilson loops. The gravity dual background is given by a particular decoupling 
limit of non-extremal (D1,D3) bound state of type IIB string theory. We obtain 
the velocity dependent quark-antiquark potential and numerically study how the 
dipole length and the potential change with velocity (for < v < 1, i.e., the Wilson 
loop is time-like) of the dipole as well as noncommutativity. We discuss and com- 
pare the results with the known commutative results. We also obtain an analytic 
expression for the screening length when the rapidity is large and the noncommu- 
tativity parameter is small with the product remaining small. When v — > 1, the 
time-like Wilson loop becomes light-like and in that case we obtain the form of 
the jet quenching parameter for the strongly coupled noncommutative Yang- Mills 
plasma which matches with our earlier results obtained using different approach. 
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1 Introduction 



One of the remarkable features of AdS/CFT correspondence [HE], [3] and its generalizations 
[I] is that it gives us access to the non-perturbative regimes of large N gauge theories 
simply from the low energy, weakly coupled string theory in certain backgrounds. So, for 
example, the expectation values of Wilson loops, which are non-perturbative objects in 
gauge theories, can be computed using AdS/CFT correspondence as has been prescribed 
in El [71 |8]. In strongly coupled gauge theories of interacting quark-gluon plasma, 
Wilson loops can be related to various measurable quantities in heavy ion experiments in 
RHIC or in LHC For example, the expectation value of a special time-like Wilson loop 
can be related to the static quark-antiquark potential (9] in a moving quark-gluon plasma. 
On the other hand, the expectation value of a particular light-like Wilson loop can be 
related, among other things, to the radiative energy loss of a parton or the jet quenching 
parameter [TO] . 

The velocity dependent quark-antiquark potential of a dipole moving with an arbitrary 
velocity through the hot quark-gluon plasma including the screening length [TT} [121 [El EH] 
as well as the jet quenching parameter [151 HBp have been calculated when the plasma is 
described by D = 4, J\f = 4, SU(iV) Yang-Mills theory using AdS/CFT correspondence^]. 
In this paper we calculate the velocity dependent quark-antiquark potential, the screening 
length and the jet quenching parameter when the plasma is described by D = 4, thermal, 
noncommutative Yang-Mills (NCYM)@ theory at large N using gauge/gravity duality. 

4 Also see [17] for a recent review. 

5 Jet quenching parameter in various other theories have been obtained in [T5]. Also the drag force on 
a moving quark have been calculated in |19j . 

6 The space-time noncommutativity is an old idea introduced first by Heisenberg and Pauli |20j in order 
to avoid infinities in quantum field theory before the renormalization was successful. It was Snyder [21] 
and then Connes [22] who took the idea seriously. Connes along with Chamscddine |23j even introduced 
noncommutative geometry as a generalization to Ricmannian geometry and obtained gauge theory as a 
companion to general relativity giving rise to a true geometric unification. In this framework parameters 
of the standard model appear as geometric invariants. Even though a consistent gauge theory can be 
formulated in noncommutative space-time, so far, there is no evidence for its existence in nature in low 
energy. However one can not rule out the possibility that its effect could be detected in some future 
experiments. The experimental lower bound on the noncommutativity scale reported in the literature 
[2"1] usually gives a very small effect and is hard to detect. So, it is desirable to look for its effect in as 
many different cases as possible. High energy heavy ion collision is one such possible case and it may be 
worth while to look whether it can provide a better window for the effect of space-time noncommutativity 
to be observed. One might wonder how would the space-time noncommutativity appear in the heavy 
ion collision in the first place? It is known that one of the mechanisms for the appearance of spatial 
noncommutativity is the presence of an intense magnetic field in the background. It has been shown 
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NCYM theory arises quite naturally in string theory [27J EEJ [29] and M-theory [30] and it 
is of interest to see how noncommutativity affects the known velocity dependent quark- 
antiquark potential, screening length as well as the jet quenching parameter of the ordinary 
super YM theory. 

The gravity background in this case is given by a particular decoupling limit [28] of 
non-extremal (Dl, D3) bound state system of type IIB string theory. (Dl, D3) bound 
state (Hll [32] contains a non-zero 5-field and it becomes asymptotically very large in 
the decoupling limit and is the source of space-space noncommutativity [27]. We use 
the string probe in this background and extremize the Nambu-Goto string world-sheet 
action in a particular static gauge and in turn obtain the expectation value of the Wilson 
loop, where the loop is the boundary of the above minimal area. We consider both the 
time-like and the light-like Wilson loops. From the time-like Wilson loop, we calculate 
the velocity dependent quark- ant iquark potential, when the background or the plasma is 
moving with a velocity v (for < v < 1) and the end-points of the fundamental string act 
as a heavy quark-antiquark pair or a dipole in the boundary gauge theory. We will take 
the background (or the plasma) to be moving along one of the brane directions (which 
is taken to be a commutative direction) relative to the dipole which lies along one of 
the noncommutative directions. There are other possibilities, but this is the simplest 
case where the noncommutative effect shows up. The quark-antiquark potential obtained 
from the minimal area of the world-sheet can only be given numerically. We first plot the 
dipole length as a function of certain constants of motion and then using this we plot the 
potential as a function of the dipole length at different velocities and noncommutativity 
parameters. The results are compared with the known commutative case [16]. We also 
give an analytic expression for the screening length when the rapidity is large and the 
noncommutativity parameter is small with the product also remaining small. Finally, for 
the completeness, we consider the case where the velocity v — > 1 and the time-like Wilson 
loop becomes light-like. In this case we recover the form of the jet quenching parameter 
as obtained before [33] using a different approach. 

This paper is organized as follows. In section 2, we compute the time-like Wilson loop 
and give our numerical results for the dipole length and the quark-antiquark potential. 
We also obtain the analytic form of the screening length in some spacial case. In section 
3, we take the velocity to unity and obtain the light-like Wilson loop. From this we obtain 
the form of the jet quenching parameter in this theory. Finally, we conclude in section 4. 

in both analytic calculations and numerical simulations [53] that such an intense magnetic field is 
indeed possible in the heavy ion collision in RHIC (or in LHC). So, it may be quite relevant to consider 
such a possibility in the present context. 
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2 Time-like Wilson Loop and Q-Q potential 



In this section we will compute the time-like Wilson loop for the hot noncommutative 
Yang-Mills theory in (3+l)-dimensions from gauge/gravity duality. The gravity dual of 
noncommutative Yang-Mills theory is given by a particular decoupling limit [2SJ of 
the non-extremal (Dl, D3) bound state of type IIB string theory. We use the funda- 
mental string as a probe in this background and compute the Nambu-Goto string world- 
sheet action. This action is then extremized to relate it to the expectation value of the 
time-like Wilson loop [16]. We will study numerically both the separation length of the 
quark- ant iquark and the potential when we vary the velocity and the noncommutativity 
parameter. We will also give an analytic expression of the screening length in some special 
limit. 

The non-extremal (Dl, D3) bound state solution of type IIB string theory is given by 
the following metric (given in the string frame), the dilaton, the NSNS 5-field and the 
RR form-fields 



ds 2 = H~\ (-f{dtf + (da; 1 ) 2 + j ((dx 2 ) 2 + (dx 3 ) 2 )^j + H* + r 2 dQ\ 

. H tana 

** = tip* B ^ = — 

Aqi = — (H^ 1 — 1) sin a coth ip, Ah23 = — cos a coth</> + T. T. (1) 

9s 9s F 

where the various functions appearing above are, 

= 1 _ r_l h = 1+ r o sinh2 1 F = 1 + r o cos2 a sinh2 1 ( 2 ) 

Here D3-branes lie along x 1 , x 2 , x 3 and Dl-branes lie along x 1 . The angle a measures 
the relative number of Dl and D3 branes by the relation cos a = N/yjN 2 + M 2 , where 
N is the number of D3-branes and M is the number of Dl-branes per unit codimension 
two surface transverse to Dl-branes [31]. Also in the above tp is the boost parameter and 
r is the radius of the horizon of the non-extremal (Dl, D3) bound state solution. is 
the dilaton and g s is the string coupling constant. Aqi and A0123 are the RR form-fields 
corresponding to Dl-brane and D3-brane respectively. T.T. denotes a term, involving 
transverse part of the brane to make the field-strength self-dual, whose explicit form is 
not required for our discussion. B23 is the NSNS form responsible for the appearance of 
noncommutativity in the decoupling limit. 

The NCYM decoupling limit is a low energy limit for which we zoom into the region 
P], 

r < r ~ r a/ sinh (p cos a <ti r a/ sinh (p (3) 
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The above limit implies that tp is a large parameter and the angle a is close to tt/2. In 
this limit we get, 

r^sinh 2 ^ H 1 h ... 
H w — — — ~ = (4) 

r 4 F cos 2 a(l + a 4 r 4 ) cos 2 a v y 

where we have defined 

h = j-7, with, a 4 = — 5 n — (5) 

1 + a 4 r 4 r 4 sinh <p cos 2 a 

From ([T]) we notice that since the asymptotic value of £>-field is tana and a — > tt/2 in the 
decoupling limit, the 5-field becomes very large. The non-vanishing component of the 
.B-field is B23 which gives rise to a magnetic field in the D3-brane world-volume and is 
responsible for making x 2 and x 3 directions noncommutative [35]. Using (j3J), we rewrite 
the metric in (TJP) as, 

ds 2 = Q ^~~r (-fdt 2 + (dx 1 ) 2 + h Udx 2 ) 2 + (dx 3 ) 2 ]) + r ° 2si f t ( d ^ + r 2 dtti) (6) 

Here the function h is as defined in flHJ) and also in writing (|6]) we have rescaled the 
coordinates as, x 2,3 — > cos ax 2,3 . This metric along with the other field configurations 
given in ([1]) in the NCYM decoupling limit is the gravity dual of four dimensional thermal 
NCYM theory. We use fundamental open string as a probe and consider its dynamics 
in this background. Let the line joining the end points of the string or the dipole lie 
along x 2 one of the noncommutative directions and move along x 1 with a velocity v where 
< v < 1 0. To simplify the calculation we can go to the rest frame (t', x 1 ') of the dipole 
by boosting the coordinate as, 

dt = cosh r]dt' — sinh rjdx 1 ' 
dx 1 = — sinh rjdt' + cosh rjdx 1 ' (7) 

where the rapidity r\ is related to v as tanhr/ = v. So, in this frame the dipole is static and 
the background (or the plasma) is moving with a velocity v in the negative x 1 -direction. 
Note that the rectangular Wilson loop lies along t' and x 2 directions and we denote the 



7 There are various other possibilities one can consider, for example, the dipole lies along the com- 
mutative direction x 1 and moves along one of the noncommutative directions x 2 (say) or the dipole lies 
along one of the noncommutative directions x 2 and moves along the other noncommutative direction x 3 . 
Dipole can even have an arbitrary orientation with respect to its motion and the motion can also be 
in arbitrary direction in the mixed commutative-noncommutativc boundary. Here we consider only the 
simplest case to see the noncommutative effect. 
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lengths along those directions as T and L respectively. We also assume T ^> L such that 
the string world-sheet is time translation invariant. Using (J7|) in the metric fl6]) we get, 

2 7 

ds 2 = -A(r)dt 2 - 2B(r)dtdx 1 + CirMdx 1 ) 2 + / , \(dx 2 ) 2 + (dx 3 ) 2 ] 

Tq smh ip 

Jl^dr 2 + 



r 2 f 



G^dx^dx 1 



where 



A{r) 
B(r) 



r ( r^cosh 77 



r\ sinh ip \ r 4 
Tq sinh 77 cosh 77 
r 2 sinh 



Note that since we will be using the 'primed' coordinates from now on, we have dropped 
the prime for simplicity. We will evaluate the Nambu-Goto action of the string world-sheet 
in this background. The Nambu-Goto action is given as, 

1 



S = J dadTy/-detg aj 3 (10) 

where g a p is the induced metric on the world-sheet given as 

Ox^ 1 0x u 

g a p = g ^q^q^ ( n ) 

with is as given in (JSJ). In the above £ a are the world-sheet coordinates r and a for 
a = and a = 1 respectively. We choose the static gauge for evaluating (ITU]) as, r = t, 
a = x 2 , where —L/2 < x 2 < L/2 and r = r(cr), x l (a) = x 3 (a) = constant. r(a) is the 
string embedding we want to determine by extremizing the Nambu-Goto action, with the 
boundary condition r(±L/2) = tqA (here A is a parameter). Using these in the action 
(ED, we get, 



s T 



2ira' 



rL/2 




/ da 


A(r) | 


/-L/2 





r 2 /i Tq sinh 99 (c^r) 



r 2 . sinh ip r 2 f 



(12) 



with -A(r) as given in (Q. Introducing dimensionless quantities y = r /r , a = a/ (ro sinh ip) 
and £ = L/(r Q sinh ip), the Nambu-Goto action flT2|) can be rewritten as, 

S=-^ daC = TTVX / (13) 

vra' ./n ./n 
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where 

y,2 



C = J (V - cosh 2 rj) [ ir^ + -r (14) 

Note that we have removed 'tilde' from a while writing (!T3|) as it is an integration variable. 
Also y' denotes dy/da and we have used the fact that y is an even function of a by 
symmetry. In writing the second expression in (fTBl we have made use of the standard 
gauge/gravity relation 



r 4 sinh 2 <p = 2g\ M Na' 2 = Ac/ 2 



rrro cosh ip titq sinh ip 

and a 4 r 4 = ^ — = vr 4 AT 4 # 2 (15) 

sinh (p cos 2 a 

The first relation in (fT5"j) is obtained from calculating the Hawking temperature of the 
non-extremal (Dl, D3) brane given by the metric in ([T|) (note that in the decoupling limit 
when ip is large cosh ip pa sinh ip) and this is the temperature of the NCYM theory by 
gauge/gravity duality. The second relation is obtained from the D3-brane charge where N 
is the number of D3-branes and in NCYM theory this is the rank of the gauge group, gyyi 



is the NCYM coupling and 2g\ u N = A 13 is the 't Hooft coupling of the NCYM theory. 
The NCYM 't Hooft coupling is related to the ordinary 't Hooft coupling by the relation 
A = (a'/9)X, where 9 is the noncommutativity parameter defined by [x 2 , x 3 ] = i9 [28] . 
Here 9 is a finite parameter and in the decoupling limit as a' -)■ 0, A remains finite. The 
third relation is obtained using the first two and also using cos a = a'/ 9 in the decoupling 
limit. Note that in the decouping limit as a' — > 0, a — > it /2 as mentioned earlier. We will 
compute y{u) by extremizing the action (fTBl . 

Now since the Lagrangian density in (JT3l does not explicitly depend on a, we have 
the following constant of motion, 

n, r fdC y A - cosh 2 rj 

ft = L — y —— = ; = q = const. (16) 

dy' 



l + aV^V^-cosh 2 ^) 



As in the commutative theory [16] we will consider two cases: (a) \Jcos\i rj < A and then 
take A — > oo. The rapidity in this case remains finite, the Wilson loop is time-like and 
the action is real. We will compute the quark-antiquark potential in this case and also 



give an expression of screening length in some special case, (b) -y/cosh r\ > A and then 
take 77—7-00, keeping A finite. The Wilson loop in this case is light-like and the action is 



3 We have used the convention adopted in JTSJ [16 
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imaginary. We will take A — > oo in the end and obtain the expression of the jet quenching 
parameter for the noncommutative hot Yang-Mills plasma. 

We will consider case (a) in this section and case (b) in the next section. When 
y^cosh 7] < A, the action would be real and from (1161) y' can be solved as, 



y' = i,\ a X°t vV - 1) (y i - vt) (it) 

q(l + a 4 r 4 y 4 ) 

where y\ = (cosh 2 1] + q 2 )/(l — a 4 r 4 g 2 ) > 1 denotes the larger turning point where y' 
vanishes. Integrating ( TTTj) we obtain, 

1/2 In r A 14- o 4 r 4 ?/ 4 

da = i{q) = — 2q / l + ar y 

y/1 - a 4 r V J yc ^^-l)^-y±) 



We remark that if we naively take A, where the boundary theory is supposed to live 
to oo, the above integral diverges. Here I is related to the dipole length L by I = 
L/(ro sinh <p) = irLT and so the divergence in £(q) is physically meaningless. Note that 
£(q) in the commutative theory is indeed finite as can be seen from f[T8"j) by putting a 2 r 2 ~ 
6, which is a measure of noncommutativity, to zero. However, for the noncommutative 
case £(q) is divergent if we take A — > oo. The reason is that the noncommutative gauge 
theory does not live at A = oo, but at some finite value of A, whose exact value is not 
known. This is implicit in [28]. It has been shown there that in a noncommutative 
theory it is not possible to fix the position of the string end point at infinity as a small 
perturbation would change it violently. This has also been explicitly mentioned in [36J. 

In the context of a Wilson loop calculation, it has been noticed before |37] that the 
string end points for a static string can not be fixed at a finite length at A — > oo in a 
noncommutative theory. Therefore, the dipole length L indeed diverges. The reason for 
this divergence has been argued to be the non-local interaction between the Q-Q pair 
in a magnetic field [38J. To be precise, the interaction point in terms of the center of 
mass coordinate gets shifted by a momentum dependent term. Thus if the only non-zero 
component of the B- field is B23, as in our case, then by placing the dipole along x 2 , it 
automatically gets a momentum along x 3 . So, if we keep the dipole static along x 3 , the 
length will diverge at infinity. To compensate the momentum along x 3 , the dipole must 
move along x 3 with a particular velocity (39]. In that case, the end points of the string 
can be fixed at a finite length on the boundary at infinity and thus the divergence in the 
dipole length gets removed. 

In the following, we, however, remove the divergence in the dipole length in a different 
way. We have noticed in (Tl8l) . the integral representation of the dipole length, that by 
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explicit evaluation the integral diverges (note that the dipole in our case does not move 
along x 3 ) as we take A — > oo. But, the integral can be regularized to give a finite result 
if we can explicitly extract the divergent term (the part that goes to infinity as A — > oo) 
of the integral and remove no After regularization, the noncommutative theory may 
be thought of as living at A = oo. So, after regularization we will take A — > oo. It 
is not difficult to obatin the divergent term from the integral (1T8|) when A — > oo and 
by inspection it can be seen to have the form 2ga 4 rgA/ a/1 — cfir\q 2 . Subtracting the 
divergent part the finite £(q) can be written in the integral form as, 



2q 



^fl~^ a*r$q 



4»,4„2 



i i 4 4 4 

1 + a 4 r 4 ?/ 4 4 4 

-ay — a r A 



lie 



V(y 4 -V(y 4 -yt) 



(19) 



The above equation therefore gives us the quark antiquark separation L(q) = £(q)/(7rT) 
of the dipole as a function of the constant of motion q. It is not possible to perform 
the integration in ( |T9l) and give an analytic expression of £{q). So, we will perform the 
integration numerically for different fixed values of the rapidity r\ and the noncommuta- 
tivity parameter ar ~ X-^TyO and use this to evaluate the quark- antiquark potential. 
However, we can give the analytic expression of £{q) and the screening length only for 
some special values of the rapidity and the noncommutativity parameter which will be 
discussed later. 

Now substituting y' from (fTTj) into the action (1T31) we get, 

s(f) = tt ^ r z 4 --" 2 " -y ( 20) 

Vi - a>ru* 4 V(y'- l )(v'-vi) 

However, as it is known [16] from the commutative case that this action is divergent as we 
take A — > oo. The reason for this divergence for the commutative case is that this action 
contains the quark- antiquark self-energy Sq and once this is subtracted from S(£), the 
action becomes finite and from there one calculates the potential. In the noncommutative 
case also, we need to subtract the quark- antiquark self energy Sq. This is calculated by 
evaluating the Nambu-Goto action for the single string stretched between the horizon tq 
and the boundary at r A and then multiply by 2 (for the two strings associated with 



9 We have checked explicitly, by doing the Wilson loop calculation in the zero temperature and zero 
rapidity (ij — > 0, implying that the dipole is not moving along a; 1 ) case and comparing it with eq.(45) 
of [37] that the two ways of regularizing (namely, by physically moving the dipole along x 3 with a 
particular velocity [39) or by removing the divergent part of the integral of £{q)) the dipole length in the 
noncommutative theory give identical results when the turning point (y c ) is large. However, for small 
turning point dipole lengths obtained by these two methods differ by some finite term. The same is true 
for the calculation of Q-Q potential. We would like to thank the anonymous referee for raising this issue. 
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quark and antiquark). The Nambu-Goto action in this case is evaluated by considering 
a quark moving in the x 1 -direction (which is commutative) and so, the gauge condition 



one uses is r = t, a = r, x 1 = x 1 (a), x 2 (a) = x 3 (a) = constant! 10 !. However, since x 2 
and x 3 are the noncommutative directions the computation of Sq remains unaffected by 
the noncommutativity. Therefore, Sq takes exactly the same form as evaluated in the 
commutative case given as [16] 



A 

A J dy 



(21) 



So, subtracting So from S(£) we get 

ttVx 



S(£) - S 



y/l — a 4 r 4 g 2 



dy 



cosh 2 rj 



— \ 1 1 — a 4 r 4 g 2 



■y/l-a*rtq*(y c -l) 



(22) 



It is clear from ( 122]) that in the commutative case when a 2 ?" 2 . ~ 9 is put to zero, S(£) — So is 
indeed finite when we take A — > oo. However, this is not the case when noncommutativity 
is present. The reason is, as we mentioned before, the NCYM theory does not live at 
A = oo, but at a finite A [281 EE]. (Another way of understanding this divergence has 
been mentioned earlier from the arguments given in [37J [381 EH]-) So, we will subtract 
the divergent term from (122!) and then take A — > oo. Doing that we find the finite 
quark- antiquark potential as, 



E(£) 



S — Sq — S { 



div 



r 



a^r^q 2 



dy 



y A — cosh r) 

vV-i)0/ 4 -£> 



1 - a%q 2 (y c ^l-a^q 2 ^ A 



(23) 



A— >oo 



where in the above Sdw is (1 — \A — a 4 r^g 2 )A with A — y oo. Here also it is not possible to 
perform the integration in (J2~3l in a closed form. So, we will obtain the quark- antiquark 
potential numerically. We first plot £(q) vs. q using (TT9T) and use it to plot E(£) vs. £ 
from (1231) at different fixed values of the rapidity 77 and the noncommutativity parameter 
ar . 



10 By constant here we mean that x 2 and x 3 are not functions of r = a. However since these are the 
noncommutative directions they will have the expected fuzziness in the x 2 -x 3 plane. 
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2.1 Plots and discussion of the results 



In this subsection we give and discuss the various plots of quark-antiquark separation 
£(q) as a function of constant of motion q and the velocity dependent quark-antiquark 
potential E(£) as a function of the quark-antiquark separation length £ for various values 
of the rapidity (77) as well as the noncommutativity parameter (ar ~ VO). 




(a) 



til 



— • 

-I 

-2 
-3 
-4 
-5 
-6 

-7 / 



ri=o.i 



ar =0.0 
ar =0.1 
ar =0.2 
ar =0.5 
ar n =1.0 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
/ 

(b) 



Figure 1: (a) shows the plot of the quark-antiquark separation l{q) as a function of the const, 
of motion q for different values of the noncommutativity parameter aro, when the rapidity is 
kept fixed at r/ = 0.1. (b) shows the plot of the properly normalized quark-antiquark potential 
E(£) as a function of £ for the same set of values of the noncommutativity parameter with the 
same value of 77 = 0.1. 



Note that each figure contains two parts (a) and (b). While in (a) we plot £(q) vs 
q, in (b) we plot the quark-antiquark potential E(£) vs £. In plot (b) we make use of 
plot (a) to obtain q for each value of £{q). Generically in all these figures the £{q) in 
part (a) starts at zero for q = (in fact for small q, £(q) — > as q which can be seen 
from ( f!9l) ) and as q increases it ends at zero either at a finite value of q, denoted as 
<7max; (when noncommutativity is present or ar 7^ 0) or at q — > 00 (when there is no 
noncommutativity or aro = 0) (for large q, £(q) — > as q~ x l 2 which can be seen again 
from ([19]) ). In between, £{q) has a single maximum at some finite q < g max indicating that 
there exists a screening length (proportional to £ max ) beyond which there is no solution to 
( TT9l) . So, given a curve, there is an £ max beyond which there is no dipole solution, that is, 
the dipole dissociates and below £ max there are two dipoles at a fixed £ for two different 
values of q. Accordingly, E{£) or Q-Q potential in part (b) generically in all the figures 
have two branches corresponding to the dipoles with two different values of q. The smaller 
value of q corresponds to the upper branch and the dipole has higher energy, whereas, 
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Figure 2: (a) shows the plot of the quark-antiquark separation £(q) as a function of the const, of 
motion q for different values (but larger than those given in Figure 1) of the noncommutativity 
parameter aro with the same r/ = 0.1 as in Figure 1. (b) shows the plot of the properly 
normalized quark-antiquark potential E{£) as a function of I for the same set of values of the 
noncommutativity parameter with the same rj = 0.1. 



the larger value of q corresponds to the lower branch and the dipole has lower energy. So, 
the dipole with lower q will be metastable and will go to the state with higher q as it is 
energetically more favorable. 

Also from part (b) we note that there exists a critical value of the rapidity r] c (whose 
value changes with the value of the noncommutativity parameter), beyond which the 
entire upper branch of the E{£) curve is negative. However, below this value, i.e., 77 < rj c , 
the E(£) curve crosses zero at £ = £ c , continues to rise till I = £ max and then turns back 
crossing zero again at I = £' c > £ c . There are various portions of the E(£) vs £ curve, 
namely, £ < £ e , £ — £ c , £ c < £ < £' & £ = £' c , £' c < £ < £ max and i > £ max which need 
separate discussions as far as the quark-antiquark potential is concerned. However, the 
discussion is pretty much the same as in the commutative case given in [161 HO] and will 
not be repeated here. We will point out the differences due to noncommutativity as we 
go along. 

In Figures 1 and 2, the rapidity is fixed at 77 = 0.1. However, in Figure 1 the non- 
commutativity parameter ar takes small values starting from 0.0 (where there is no 
noncommutativity) to 1.0, whereas, in Figure 2 the noncommutativity parameter takes 
fairly large values starting from 2.0 to 10.0. The main difference between the commu- 
tative results and the noncommutative results is that in the former case the constant of 
motion q can take arbitrarily large values, but in the latter case q can not exceed certain 
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Figure 3: (a) shows the plot of the quark-antiquark separation l(q) as a function of the const, 
of motion q for different values of the noncommutativity parameter aro, but now the rapidity is 
kept fixed at r/ = 0.5. (b) shows the plot of the properly normalized quark-antiquark potential 
E(£) as a function of £ for the same set of values of the noncommutativity parameter with the 
same value of r\ = 0.5 



finite value (g max ) because beyond this value the quark-antiquark separation £(q) becomes 
negative which is unphysical. The reason behind this cut-off is the regularization of the 
integral made in ( !T9l) . Note that the last term in ( !T9|) is subtracted in order to make 
£(q) finite as A — > oo. However, as q increases, y c increases which makes the last term 
dominate over the integral and therefore £(q) becomes negative. Thus this effect is due 
to the noncommutativity of the underlying boundary theory. We see from Figure 1(a) 
that as the noncommutativity parameter aro increases, £(q) curve deviates more and more 
from the commutative curve, the £ max falls and the peak shifts towards the left (i.e., the 
maximum occurs at a smaller value of q). In particular, the deviation from the commuta- 
tive case becomes more pronounced after £ max is reached. However this feature continues 
upto certain value of ar ~ 2.0 and as it is increased further (see Figure 2(a)) the £(q) 
curve now deviates more from the commutative case throughout the allowed range of q, 
but the maximum value, £ max , again starts rising and the peak as before shifts further 
towards left, i.e., towards smaller values of q. Note that the scales in Figures 1 and 2 are 
quite different and are chosen so such that we can see the plots for various higher values 
of ar . Figure 1(b) and 2(b) shows the plot of the velocity dependent Q-Q potential 
E(£) vs the quark-antiquark separation length £ for t] = 0.1 with various values of the 
noncommutativity parameter ar . As we mentioned before each curve in this case has two 
branches corresponding to the two dipole solutions obtained in Figures (a). The slight 
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Figure 4: (a) shows the plot of the quark-antiquark separation l(q) as a function of the const, 
of motion q for different values of the noncommutativity parameter aro, but now the rapidity is 
fixed at r/ = 1.0. (b) shows the plot of the properly normalized quark-antiquark potential E(£) 
as a function of £ for the same set of values of the noncommutative parameter with the same 

T) = 1.0. 



deviation of £(q) from the commutative case for small values of q, i.e., below the value of 
q corresponding to £(q) = £ max , (see Figure 1(a)) is reflected in that fact that in Figure 
1(b) the upper branches almost merge with the commutative counterpart whereas the 
greater deviation in £(q) after £ max is reached leads to a rise in the lower branch of the 
E{£) curve from the commutative case in Figure 1(b). However, as the noncommutativity 
parameter is increased the overall deviation (particularly in the lower branch) of the E(£) 
curve is more pronounced from its commutative value. In contrast, in Figure 2(b) as the 
noncommutativity parameter is further increased, E(£), in general, dips slightly for both 
the branches. 

The feature that the screening length (~ £ mauX ) initially drops and then rises and cor- 
respondingly the lower branch of the potential E(£) rises and then drops as we go on 
increasing the noncommutativity parameter aro (with the transition occurring at around 
aro — 2.0), occurs only for the smaller value of the rapidity, rj = 0.1. As the rapidity 
becomes higher its effect starts to dominate and the transition (from falling £ max to ris- 
ing £ max as the noncommutativity parameter is increased) is suppressed so that now the 
screening length continuously drops and the lower branch of the Q-Q potential contin- 
uously rises. We have seen this to happen for rj = 0.5 and rj = 1.0. That is why we 
have given those plots only for the smaller values of ar in Figures 3 and 4. Although 
the details of these plots are different, the general features remain very similar to those 
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q I 

(a) (b) 

Figure 5: (a) shows the plot of quark-antiquark separation £(q) as a function of the const, of 
motion q for different values of the rapidity 77 with the noncommutativity parameter kept fixed 
at aro = 0.1. (b) shows the plot of the properly normalized quark-antiquark potential E{£) as 
a function of £ for the same set of values of the rapidity with the same aro = 0.1. 



of 77 = 0.1 as we have discussed and so, we do not elaborate their discussion further to 
avoid repetitions. 

We have pointed out that with different rapidities the general features of the plots 
1, 3 and 4 (i.e., for small noncommutativity) are quite similar although the details are 
different. We see from plots (a) that in all three cases, the maximum of the £(q) curve 
(l max ) or the screening length drops as the noncommutativity is increased. This implies 
that with the increase of noncommutativity, less and less number of dipoles will be formed 
in the QGP and there will be more J /\& suppression [H]. On the other hand, from plots 
(b) we see that with the increase of noncommutativity the Q-Q potential rises (see the 
lower curves which correspond to the stable states) in value which means that the quark 
and antiquark will be more and more loosely bound and eventually there will be no bound 
state formation. This may be expected since there is a fuzziness in the direction of the 
dipole due to noncommutativity. However, in plot 2(a), when the noncommutativity is 
large (and the rapidity remains small) we see that around ar = 2.0, the maximum of the 
£(q) curve (£ ma x) starts rising again and so more dipoles can form, but from plot (b) we see 
(from the lower curve) that in this case the quarks and antiquarks will be very very loosely 
bound. This does not happen when the rapidity is large (we have not shown the plots for 
this case with large noncommutativity). Here the screening length (~ ^ max ) continually 
drops as we have checked. This seems to indicate that the effect of noncommutativity 
(for large noncommutativity) gets suppressed when the rapidity or the velocity of the 
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Figure 6: (a) shows the plot of the quark-antiquark separation l(q) as a function of the const, 
of motion q for different values of the rapidity rj but now the noncommutativity parameter is 
kept fixed at a much larger value aro = 10.0. (b) shows the plot of the properly normalized 
quark-antiquark potential E{£) as a function of t for the same set of values of the rapidity with 
the same aro = 10.0. 



dipole is high. But for small noncommutativity the effect is quite prominent even when 
the velocity of the dipole is high (as shown in plots 3, 4). 

In contrast to Figures 1-4, where we plot £(q) vs q and E(£) vs £ for fixed value 
of the rapidity rj but with varying values of noncommutativity parameter aro, in Figures 

5, 6, we plot the same functions for fixed value of noncommutativity parameter ar , but 
varying values of rapidity rj. In Figure 5, ar is fixed to a small value 0.1 whereas in Figure 

6, it is fixed to a large value 10.0. In Figure 5(a) we find that as the rapidity increases 
the screening length decreases (which means there will be less dipole formation i.e., more 
J '/ty suppression) and the peaks shift towards right, i.e., to the larger value of q. This is 
expected as in the commutative case also there is a decrease in screening length with the 
increase in rapidity. Further note that for larger value of q, £(q) becomes independent 
of q. This is also shown in E{£) vs I plot in Figure 5(b), i.e. the lower branches of the 
curves merge. Contrast this to the case when aro is changed (but still kept small) keeping 
rj fixed, when the lower part of the £(q) curve (i.e., small q) does not exhibit significant 
deviation, the peak shifts towards left and the upper branch of the E(£) curves merge. 
So we can think of rj and aro as sort of having opposite effects. These features are also 
evident for large values of aro given in Figure 6. In figure 6(a) as the rapidity rj increases 
the screening length decreases and the peaks shift towards right, but since now the scale 
of g-axis is very much enlarged this is not much visible (this is also due to the fact that rj 
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now changes by a very small amount). The independence of £(q) with q for larger values 
of q is not evident in this case due to the differences in scale in g-axis for Figure 5(a) and 
Figure 6(a), however, it is clear that it has this tendency. Unlike in Figure 5(b) the lower 
branches of the E(£) curves do not merge for different values of r\ as is evident from Figure 
6(b). However, the spread is again due to the enlarged (compared to Figure 5(b)) scale 
of the E(£) axis which is chosen to show the two branches of the E(£) curve distinctly. 



2.2 Screening length in a special case 

The expression for the regularized quark-antiquark separation length £(q) as a function of 
the constant of motion q is given in f fT9|) . However, as we have mentioned, it is not possible 
to perform the integration occurring there in a closed form in general and give an exact 
analytic expression for £(q). That is why we have numerically solved ( f!9|) and plotted 
£(q) in the previous subsection. This has nothing to do with the noncommutativity of the 
underlying gauge theory and this happens also for the case of commutative theory. For 
the case of commutative theory it is possible to give an exact analytic expression of £(q) 
only in the large velocity or large rapidity limit. Noncommutativity, on the other hand, 
makes the analysis a little bit more involved and in this case it is possible to obtain the 
analytic expression only when the rapidity is large and the noncommutativity is small 
with the product remaining small. For large rj or large y c , the expression for £(q) in f fT9|) 
can be expanded as follows, 



£{q) 



2q 



4 4 4 



y/l - a^q 2 J yc y 2 Jt^Vi 
3q 



dy + 



i i 4 4 4 

1 + a r^y 



y/1 - a^q 2 J yc y«y/y*-y. 



1 + a 4 r 4 y 4 2ga 4 r 4 



Vl - a%y J yc y^^/y^^t 



dy 



dy 



:A 



'1 - a 4 r 4 g 2 



(24) 



When A — > oo, the above integrals can be evaluated and £(q) can be written as a series 
expansion in inverse powers of y c as, 



£{q) 



1 — a^r^q 2 



-2a\ 4 ^ + (2 + aS 



r (!) vt 



4r (1) yl 



(25) 

By construction the divergent last term in eq. (T24"]) gets cancelled with the divergent term 
in the first integral when A — y oo. The other integrals are convergent and makes the 
expression for £(q) finite. By taking the first three terms in the series we can obtain the 
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value of g and y c which maximizes £(q) as, 



2 cosh r/(l — 15a r cosh 77) 



cosh r] + q 2 
(1 - a 4 r 4 g 2 ) 



3 cosh 77(1 — 8a r cosh 77) 



(26) 



In obtaining the above expressions we have assumed a 4 r 4 <C 1 and a 4 r 4 cosh 2 77 <C 1. 
Using ( 1261) we obtain the maximum value of £ upto next to leading order as, 



2^ (f) 

3 3 / 4 r (i) COshS 77 



7 
2' 



1 — -a 4 r 4 cosh 2 77 + 



0.74333 



7 



1 — 2 a4r o cos h 2 ?7 



cosh 2 77 

By using (115j) we can rewrite £ max in ( 1271) in terms of the gauge theory parameters as 



(27) 



0.74333(1 -v 2 ) 1 * 



7 7r 4 XT A e 2 

1_ 2 1-7J 2 + 



(28) 



where we have used cosh 77 = 7 = 1/ y/Y— v 2 , with 7j being the velocity of the dipole. In 
(12 8p the term outside the square bracket is the commutative result (when we put 6 = 0) 
and represents the usual J/\l/ suppression of the high velocity quark antiquark produced 
in the QGP in the heavy ion collision observed in RHIC [TTj l4Tj . However, we note that 
the noncommutativity reduces this result due to the second term in the square bracket 
in (I2"8"j) . The quantity L max = £ max / (vrT) can be thought of as the screening length of the 
dipole since this is the maximum value of L beyond which we have two dissociated quark 
and antiquark or two disjoint world-sheet for which E(£) = 0. As the screening length 
gets smaller less and less dipoles will be created and there will be more suppression of 
quark-antiquark bound states like J/\l/. Noncommutativity makes the interaction between 
the quark and antiquark weaker due to nonlocality and that is the reason it makes the 
screening length shorter. Note that the velocity of the dipole has an opposite effect in 
the correction term due to noncommutativity, i.e., higher the velocity, lower would be the 
correction term due to noncommutativity. Also the correction term is more pronounced at 
higher temperature. We would also like to remark that noncommutativity gives a range for 
the temperature. We mentioned that the above result (128]) is valid when a 4 r 4 cosh 2 77 1 
which in turn gives a range for the temperature as, 



T < 



tt 4 A(1 -v 2 )9 2 



(29) 
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When the temperature is above this value the expansion in (127)1 will break down and the 
screening length will no longer be given by ( )28|) . In that case the screening length has 
to be computed in the opposite limit where a 4 r 4 cosh 2 rj 3> 1. However, in this limit we 
haven't been able to write a closed form analytic expression for the screening length. 



3 Jet quenching parameter 

In section 2, we have discussed the case (a) where the rapidity rj remains finite and 
-^/cosh r] < A. So, the velocity of the background is in the range < v < 1 and the 
Wilson loop is time-like. Now we will discuss the case (b) where -\/cosh r] > A. In order 
to calculate the jet quenching parameter we take r\ — > oo or v — > 1, so that the Wilson 
loop is light-like and then take A — > oo. The jet quenching parameter in the NCYM 
theory has been calculated before [33] in a different approach, namely, we calculated it 
directly in the light-cone coordinates. Here we compute it, primarily for completeness, 
from the general velocity dependent string world sheet action and then taking the limit 
v — > 1. As before we compute the expectation value of the Wilson loop by extremizing 
the action. The jet quenching parameter, a measure of radiative parton energy loss in a 
quark-gluon plasma medium, is related to a particular light-like Wilson loop and thus we 
obtain its value using the holographic gauge/gravity duality. Since this has already been 
calculated in [33], using a different approach, we will be brief here. Note that as cosh 2 1] is 
now greater than A 4 , where A is the upper limit of y, the factor (y 4 — cosh 2 77) appearing 
in the action f|T3|) . f|T^|) is negative and the action becomes imaginary. So, we rewrite the 
action ffT3l) as, 



S = — - / daC = iTTVX / daC (30) 



where 



£=J(cosh 2 ,-, 4 ) 



y,2 



a 4 r V ' y*-lj (31) 
As before since the Lagrangian density does not explicitly depend on a, the corresponding 
Hamiltonian is conserved. Therefore, we have 

v r > dC + cosh 2 rj -y 4 

ft = C - y — = const. = q (32) 

V (1 + a 4 r V) ^(cosh 2 ^ - y 4 ) + 

where we have denoted the constant as go- fl3"2"j) can be solved for y' as 



, _ y/l + a^goVd/ 4 -!)^-!/ 4 ) 
V ~ go(l + « 4 r V) {66) 
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where 

On integration, fl33|) gives us, 

rl/2 

da 



cosh 2 r] — g 2 



a% 4 g 2 



2g 



4 4 4 

1 + a .r^y 



-.dy 



'1 + a^ql A y/tf-lKfa-yi) 
Substituting the value of y' from fl33|) into the action f|3T|) . we simplify its form as, 



(34) 



(35) 



(36) 



iTTV X f A cosh 2 rj — y 4 
y/1 + a 4 r 4 g 2 A ^ _ x)^ _ y 4)' 

Now in the above since £ is related to the dipole length as L = £/(ttT), it is very small 
compared to the other length scale in the theory and so, from fl35l) it is clear that go is 
also a small parameter. In this approximation go can be obtained from fl35l) as 



go 



' cosh rj 



4 4 4 

1 + a r$y 



l 



■dy 



In this limit S(£) in ( 1361) can be expanded as, 

S(£) = S<® + q 2 SV + O (q 4 ) 



(37) 



(a 



where 



5 (0) 



2(S (i) 



-dy 



-% 2 



4 4 4 



(39) 



'i \/ cosh 2 77 — y 4 \/y 4 — 1 

It can be shown [T5], [16] that as go — > 0, above is equal to So, the self-energy of the 
dissociated quark and antiquark or area of the two disjoint world-sheet. So, subtracting 
the self-energy we obtain the action as 

-1 -1 



S - S = q 2 S^ = l H^f C osh 77 



444 
- a r^y 



dy 



(40) 



where we have used (1371) and have taken 77—7-00. Now T cosh 77 in (T40I) can be identified 
as L~/\/2, where L~ is the length of the Wilson loop in the light-like direction. Let us 
use the standard relation (for example see, [TO] ) 



(W(C)) = e 2i(5(c) " 5o) 



(41) 
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where the factor 2 in the exponent in the second expression is due to the fact that we are 
dealing with adjoint Wilson loop. Here L and ^ncym is the jet quenching parameter 
of the NCYM theory. Now we can extract its value from (14011 as, 



<?NCYM = VT 2 V / It 3 



1 i 4 4 4 

1 + a 4 r 4 ?/ 4 

-dy 



y A 



(42) 



We would like to point out that the above expression of the jet quenching parameter is 
actually a formal expression since by taking A — > oo, the integral in the square bracket in 
the jet quenching expression diverges. The reason for this divergence, as we mentioned, is 
that for the noncommutative case the gauge theory does not live at A = oo, but at some 
finite A [281 EH] (see also [37J ES [39]). So, the above integral needs to be regularized. 
We have given the details of the regularization in [33] and here we just give the results, 
namely, the regularized integral has the value, 

Substituting ( )43|) in ( )42|) and expressing a 4 r 4 in terms of the gauge theory parameters 
from ( Tl~5]) we obtain, 



3 

7T2 



qncym = ~Yjiy VXT ( 1 + — 3 — 

So, for small noncommutativity, 6 -C 1, the jet quenching parameter is given as, 



3 

7T2 



7T 4 T 4 A# 2 



<?ncym = ~~^jiy~ 1 1 §— + ) ( 45 ) 

Whereas, for large noncommutativity, 9 ^> 1, the jet quenching parameter takes the form, 

*™-3^( i -sb +o £ ) ) (46) 

We note from fT45|) that for small noncommutativity, when 9 — > 0, we recover the value 
of the jet quenching parameter of the ordinary Yang-Mills plasma obtained in [15] as 
expected. In this case the NCYM 't Hooft coupling A reduces to ordinary 't Hooft 
coupling. But in the presence of noncommutativity, the value of the jet quenching gets 
reduced from its commutative value and the reduction gets enhanced with temperature 
as T 7 . The reduction in the jet quenching for the noncommutative case can be intuitively 
understood as follows. The noncommutativity introduces a non-locality in space due to 
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space uncertainty and there is no point-like ineraction among the partons and therefore, 
the parton energy loss would be less. Also since for small noncommutativity we have 
a 4 ?"o <C 1, this gives a range for the temperature due to noncommutativity as, 



when the temperature is above this value, the jet quenching expression will no longer be 
given by In that case we have to use the expression ( jlBJ) which is valid when the 

temperature is given by the limit 



In this case the jet quenching varies inversely with temperature. 

In [33] the jet quenching for the NCYM theory has been calculated by evaluating 
the expectation value of the Wilson loop directly in the light-cone frame. But here we 
calculate the same for finite value of the rapidity rj and then take the limit i] — > oo or the 
velocity of the background v — > 1. Obviously, both methods give the same results. 

4 Conclusion 

To summarize, in this paper we have computed the expectation value of both the time- 
like and the light-like Wilson loop of the NCYM theory in (3+1) dimensions using the 
gauge/gravity duality and the Maldacena prescription. The gravity dual background for 
the NCYM theory is given by a particular decoupling limit of (Dl, D3) brane bound state 
system of type IIB string theory. The noncommutative directions were taken as x 2 and 
x 3 , two of the world- volume directions of the D3-branes. We introduced a fundamental 
string in this background as probe whose end points or the dipole consisting of a quark 
and an antiquark lie along one of the noncommutative directions, namely, x 2 -direction. 
The background or the QGP medium is assumed to move along x 1 -direction with a finite 
velocity < v < 1. This particular configuration is taken for simplicity. The expectation 
value of the Wilson loop is calculated by extremizing the world-sheet area of the F-string 
in the background whose boundary is the loop in question. 

Initially we took the velocity to be less than 1, the case in which the action is real and 
the Wilson loop is time-like. In this case we computed the quark-antiquark separation 
length I as a function of constant of motion q and using this we also computed the quark- 
antiquark potential E as a function of £. We found that unlike in commutative case, £(q) 




(47) 




(48) 
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diverges as we take the boundary theory to be living at A = oo and so we needed to 
regularize an integral to make the result finite. After this regularization we found that 
the constant of motion q can not take arbitrarily large values as in the commutative case, 
but it must have a cut-off. However, we could not get a closed form analytic expression 
for £(q) in general and therefore, obtained it numerically and plotted the function £(q) 
vs q with fixed value of the noncommutativity parameter and different values of rapidity 
and also with fixed value of rapidity and different values of noncommutativity parameter. 
We discussed the various cases and pointed out the differences with the commutative 
results. Similarly, the quark- ant iquark potential E(£) was also found to be divergent 
even after subtracting the self-energies of the quark and the antiquark unlike in the 
commutative case and a regularization was needed. Again in this case we could not 
get an analytic expression for E(£) and we have plotted the function E(£) vs £ with 
fixed value of noncommutativity parameter and varying rapidity as well as fixed value 
of rapidity and varying noncomutativity parameter. Here also we discussed the results 
and compared them with the commutative results. Then we gave an analytic expression 
for the screening length in a special limit, namely, when the rapidity is large and the 
noncommutativity parameter is small with the product remaining small and discussed 
our results. 

Finally, we discussed the case when the rapidity goes to infinity or the velocity of 
the medium approaches 1. In this case the action becomes imaginary and the Wilson 
loop becomes light-like. Using this Wilson loop we obtained the expression of the jet 
quenching parameter in the NCYM theory which was obtained before [33] using differ- 
ent approach. We obtained the expressions for the jet quenching parameter for small 
and large noncommutativity. For small noncommutativity the jet quenching got reduced 
from its commutative value by an amount proportional to X 3 ^ 2 T 7 9 2 , whereas for large 
noncommutativity the leading order term was found to be proportional to l/{\/\T6 2 ). 

Acknowledgements 

We would like to thank Munshi G Mustafa for several discussions. We would also like to 
thank the anonymous referee for raising some points, the clarifications of which, we hope, 
has improved the manuscript. 



23 



References 



[1] J. M. Maldacena, "The Large N limit of superconformal field theories and supergrav- 
ity," Adv. Theor. Math. Phys. 2, 231-252 (1998). |hep-th/971120~0] . 

[2] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 253- 
291 (1998). |hep-th/9802150] . 

[3] S. S. Gubser, I. R. Klebanov, A. M. Polyakov, "Gauge theory correlators from non- 
critical string theory," Phys. Lett. B428, 105-114 (1998). |hep-th/9802~l~0~9~] . 

[4] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri, Y. Oz, "Large N field theories, 
string theory and gravity," Phys. Rept. 323, 183-386 (2000). |hep-th/9905lTT] . 

[5] J. M. Maldacena, "Wilson loops in large N field theories," Phys. Rev. Lett. 80, 
4859-4862 (1998). |hep-th/9803002] . 

[6] S. -J. Rey, J. -T. Yee, "Macroscopic strings as heavy quarks in large N gauge 
theory and anti-de Sitter supergravity," Eur. Phys. J. C22, 379-394 (2001). 
|hep-th/9803001] . 

[7] S. -J. Rey, S. Theisen, J. -T. Yee, "Wilson-Polyakov loop at finite temperature in 
large N gauge theory and anti-de Sitter supergravity," Nucl. Phys. B527, 171-186 
(1998). lhep-th/9803135] . 

[8] A. Brandhuber, N. Itzhaki, J. Sonnenschein, S. Yankielowicz, "Wilson loops 
in the large N limit at finite temperature," Phys. Lett. B434, 36-40 (1998). 
|hep-th/9803137] . 

[9] K. G. Wilson, "Confinement of Quarks," Phys. Rev. D10, 2445-2459 (1974). 

[10] A. Kovner, U. A. Wiedemann, "Gluon radiation and parton energy loss," In *Hwa, 
R.C. (ed.) et al: Quark gluon plasma* 192-248. |hep-ph/030415l] . 

[11] H. Liu, K. Rajagopal, U. A. Wiedemann, "An AdS/CFT Calculation of Screening in 
a Hot Wind," Phys. Rev. Lett. 98, 182301 (2007). [ hep-ph/0607062] . 

[12] E. Caceres, M. Natsuume, T. Okamura, "Screening length in plasma winds," JHEP 



0610, 011 (2006). |hep-th/0 607233|. 

[13] M. Chernicoff, J. A. Garcia, A. Guijosa, "The Energy of a Moving Quark-Antiquark 
Pair in an N=4 SYM Plasma," JHEP 0609, 068 (2006). |hep-th/0607089~] . 



21 



[14] S. D. Avramis, K. Sfetsos, D. Zoakos, "On the velocity and chemical-potential de- 
pendence of the heavy-quark interaction in N=4 SYM plasmas," Phys. Rev. D75, 
025009 (2007). [hep-th/0609079| . 

[15] H. Liu, K. Rajagopal, U. A. Wiedemann, "Calculating the jet quenching parameter 
from AdS/CFT," Phys. Rev. Lett. 97, 182301 (2006). [hep-ph/0605178] . 

[16] H. Liu, K. Rajagopal, U. A. Wiedemann, "Wilson loops in heavy ion collisions and 
their calculation in AdS/CFT," JHEP 0703, 066 (2007). |hep-ph/0612168] . 

[17] J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal, U. A. Wiedemann, 
"Gauge/String Duality, Hot QCD and Heavy Ion Collisions," |arXiv: 1101 .0618 [hep- 

th]]. 

[18] A. Buchel, "On jet quenching parameters in strongly coupled non-conformal gauge 
theories," Phys. Rev. D74, 046006 (2006). |hep-th/0605178] ; E. Caceres, A. Gui- 
josa, "On Drag Forces and Jet Quenching in Strongly Coupled Plasmas," JHEP 
0612, 068 (2006). |hep-th/06061~34] ; F. -L. Lin, T. Matsuo, "Jet Quenching Pa- 
rameter in Medium with Chemical Potential from AdS/CFT," Phys. Lett. B641, 
45-49 (2006). |hep-th/0606136 |; S. D. Avramis, K. Sfetsos, "Supergravity and the 
jet quenching parameter in the presence of R-charge densities," JHEP 0701, 065 



(2007). [ hep-th/060 6190|; N. Armesto, J. D. Edelstein, J. Mas, "Jet quenching 
at finite 't Hooft coupling and chemical potential from AdS/CFT," JHEP 0609, 
039 (2006). |hep-ph/0606245] ; E. Nakano, S. Teraguchi, W. -Y. Wen, "Drag force, 
jet quenching, and AdS/QCD," Phys. Rev. D75, 085016 (2007). |hep-ph/0608274] ; 
G. Bertoldi, F. Bigazzi, A. L. Cotrone, J. D. Edelstein, "Holography and unquenched 
quark-gluon plasmas," Phys. Rev. D76, 065007 (2007). [hep-th/0702225] ; F. Bigazzi, 
A. L. Cotrone, J. Mas, A. Paredes, A. V. Ramallo, J. Tarrio, "D3-D7 Quark-Gluon 
Plasmas," JHEP 0911, 117 (2009). [ arXiv: 0909 .28651 [hep-th]]. 

[19] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz, L. G. Yaffe, "Energy loss of a heavy 
quark moving through N=4 supersymmetric Yang-Mills plasma," JHEP 0607, 013 
(2006). | hep-th/0605158| ; S. S. Gubser, "Drag force in AdS/CFT," Phys. Rev. D74, 
126005 (2006). |hep-th/0605182] ; E. Caceres, A. Guijosa, "Drag force in charged 
N=4 SYM plasma," JHEP 0611, 077 (2006). |[hep-th/0605235| ; J. Casalderrey- 
Solana, D. Teaney, "Heavy quark diffusion in strongly coupled N=4 Yang-Mills," 
Phys. Rev. D74, 085012 (2006). [hep-ph/0605~l~9~9~| ; T. Matsuo, D. Tomino, W. - 
Y. Wen, "Drag force in SYM plasma with B field from AdS/CFT," JHEP 0610, 



25 



055 (2006). |hep-th/060 7178| ; S. Roy, "Holography and drag force in thermal plasma 
of non-commutative Yang-Mills theories in diverse dimensions," Phys. Lett. B682, 
93-97 (2009). |arXiv:0907.0333l [hep-th]]: K. L. Panigrahi, S. Roy, "Drag force in a 
hot non-relativistic, non-commutative Yang-Mills plasma," JHEP 1004, 003 (2010). 
jarXiv: 1001 .29041 [hep-th]]. 

[20] W. Heisenberg and W. Pauli, "On Quantum Field Theory, (in German)," Z. Phys. 56, 
1 (1929); W. Heisenberg and W. Pauli, "On Quantum Field Theory. 2. (in German)," 
Z. Phys. 59, 168 (1930). 

[21] H. S. Snyder, "Quantized space-time," Phys. Rev. 71, 38 (1947); "The Electromag- 
netic Field in Quantized Space-Time," Phys. Rev. 72, 68 (1947). 

[22] A. Connes, in *Cargese 1987, Proceedings, Nonperturbative Quantum Field Theory* 
33-69; "Noncommutative Geometry," Academic Press, San Diego, CA (1994). 

[23] A. H. Chamseddine and A. Connes, "Universal formula for noncommutative geometry 
actions: Unification of gravity and the standard model," Phys. Rev. Lett. 77, 4868 
(1996). 

[24] S. M. Carroll, J. A. Harvey, V. A. Kostelecky, C. D. Lane and T. Okamoto, "Non- 
commutative field theory and Lorentz violation," Phys. Rev. Lett. 87, 141601 (2001) 
|hep-th/0105082] . 

[25] D. E. Kharzeev, L. D. McLerran and H. J. Warringa, "The Effects of topological 
charge change in heavy ion collisions: 'Event by event P and CP violation'," Nucl. 



Phys. A 803, 227 (2008) |arXiv:0711.0950l [hep-ph]]. 

[26] V. Skokov, A. Y. .Illarionov and V. Toneev, "Estimate of the magnetic field strength 
in heavy-ion collisions," Int. J. Mod. Phys. A 24, 5925 (2009) jarXiv: 0907. 13961 [nucl- 

th]]. 

[27] N. Seiberg and E. Witten, "String theory and noncommutative geometry," JHEP 
9909, 032 (1999) |hep-th/9908142] . 

[28] J. M. Maldacena and J. G. Russo, "Large N limit of noncommutative gauge theories," 
JHEP 9909, 025 (1999) |hep-th/9908134] . 

A. Hashimoto and N. Itzhaki, "Noncommutative Yang-Mills and the AdS / CFT 
correspondence," Phys. Lett. B 465, 142 (1999) [ hep-th/9907166] . 



26 



[30] A. Connes, M. R. Douglas and A. S. Schwarz, "Noncommutative geometry and matrix 
theory: Compactification on tori," JHEP 9802, 003 (1998) [hep-th/9711162] . 

[31] J. C. Breckenridge, G. Michaud and R. C. Myers, "More D-brane bound states," 
Phys. Rev. D 55, 6438 (1997) |hep-th/961ll74| |; J. G. Russo and A. A. Tseytlin, 
"Waves, boosted branes and BPS states in m theory," Nucl. Phys. B 490, 121 (1997) 
hep-th/96110~47| ; M. S. Costa and G. Papadopoulos, Nucl. Phys. B 510, 217 (1998) 



|hep-th/9612204| . 

[32] R. -G. Cai and N. Ohta, "Noncommutative and ordinary super Yang- Mills on (D(p - 
2), D p) bound states," JHEP 0003, 009 (2000) |hep-th/00012l3] . 

[33] S. Chakraborty and S. Roy, "Calculating the jet quenching parameter in the plasma 
of NCYM theory from gauge/gravity duality," Phys. Rev. D 85, 046006 (2012) 
|arXiv:1105.3384l [hep-th]]. 

[34] J. X. Lu and S. Roy, "((F, Dl), D3) bound state and its T dual daughters," JHEP 
0001, 034 (2000) [hep-th/99050l4] ; J. X. Lu and S. Roy, "(p + l)-dimensional 
noncommutative Yang-Mills and D(p - 2)-branes," Nucl. Phys. B 579, 229 (2000) 
|hep-th/991216 5|. 

[35] F. Ardalan, H. Arfaei and M. M. Sheikh- Jabbari, "Dirac quantization of open strings 



and noncommutativity in branes," Nucl. Phys. B 576, 578 (2000) [hep-th/9906161] ; 
C. -S. Chu and P. -M. Ho, "Constrained quantization of open string in background B 
field and noncommutative D-brane," Nucl. Phys. B 568, 447 (2000) [hep-th/9906~l~9~2] . 

[36] A. Dhar and Y. Kitazawa, "Wilson loops in strongly coupled noncommutative gauge 
theories," Phys. Rev. D 63, 125005 (2001) [ hep-th/0010256] . 

[37] M. Alishahiha, Y. Oz and M. M. Sheikh- Jabbari, "Supergravity and large N non- 
commutative field theories," JHEP 9911, 007 (1999) |hep-th/990921~5] . 

[38] D. Bigatti and L. Susskind, "Magnetic fields, branes and noncommutative geometry," 
Phys. Rev. D 62, 066004 (2000) [hep-th/9908056] . 

[39] S. S. Haque and A. Hashimoto, "Mass-spin relation for quark anti-quark bound 
states in non-commutative Yang-Mills theory," Nucl. Phys. B 829, 555 (2010) 
|arXiv:0903.484H [hep-th]]. 



27 



[40] S. Chakraborty and S. Roy, "Wilson loops in (p+l)-dimensional Yang-Mills the- 
ories using gravity/gauge theory correspondence," Nucl. Phys. B 850, 463 (2011) 
|arXiv:1103.1248l [hep-th]]. 

[41] S. Digal, P. Petreczky and H. Satz, "String breaking and quarkonium dissociation 
at finite temperatures," Phys. Lett. B 514, 57 (2001) |hep-ph/0105234| ; F. Karsch, 
D. Kharzeev and H. Satz, "Sequential charmonium dissociation," Phys. Lett. B 637, 
75 (2006) |hep-ph/0512 239|; H. Satz, "Quarkonium Binding and Dissociation: The 
Spectral Analysis of the QGP," Nucl. Phys. A 783, 249 (2007) |hep-ph/0609T97| . 



28 



